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ABSTRACT 

We report the first science observations and results obtained with the "extended" SMA (eSMA), 
which is composed of the SMA (SubmiUimeter Array), JCMT (James Clerk Maxwell Telescope) and 
CSO (Caltech SubmiUimeter Observatory). Redshifted absorptions at z=0.886 of C I i^Pi —^Pq) were 
observed with the eSMA with an angular resolution of 0'.'55x0'.'22 at 1.1 mm toward the southwestern 
image of the remarkable lensed quasar PKS 1830—211, but not toward the northeastern component 
at a separation of 1". Additionally, SMA observations of CO, ^^CO and C^^O (aU J=4-3) were 
obtained toward this object: CO was also detected toward the SW component, but none of the 
isotopologues were. This is the first time [C I] is detected in this object, allowing the first direct 
determination of relative abundances of neutral atomic carbon to CO in the molecular clouds of a 
spiral galaxy at z > 0.1. The [C I] and CO profiles can be decomposed into two and three velocity 
components respectively. We derive C/CO column density ratios ranging from <0.5 (representative 
of dense cores) to ^2.5 (close to translucent clouds values). This could indicate that we are seeing 
environments with different physical conditions or that we are witnessing chemical evolution of regions 
where C has not completely been converted into CO. 

Subject headings: quasars: individual (PKS 1830) — quasars: absorption lines — galaxies: ISM — 
ISM: abundances 



1. INTRODUCTION 

A powerful way to study the physical and chemical con- 
ditions in the interstellar medium (ISM) of distant galax- 
ies is to observe chemical species in absorption against 
a strong continuum background source. Few such lines 
of sight have been unveiled so far but PKS 1830—211 
(hereafter PKS 1830) is one of the most remarkable sys- 
tems allowing this kind of study. PKS 1 830 is a radio- 
loud quasar with a redshift of z = 2.507 ijLidman et al.l 
|1999| ). for which millimeter continuum emission images 
show two compact components, nor th-east (NE) ari d 
south-west (SW), separated by 1" (|Frve et al.lll997D . 
This double structure represents two images of the back- 
ground quasar, magni fied and distorted by a lens ing sys- 
tem at z = 0.88582 (jWiklind fc Combed [19961 ). This 
lens was at first indirectly detected by the observations 
of broad H I and molecular absorptions at millimeter 
wave l engths dW iklind fc Combes 1996, 1998:' Gerin et all 
[19971 iMuUer et al... ,2006.1 . Wiklind fc Combed (119981 ) 
suggested the lensing system to be an almost face-on 
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spiral galaxy, which was recently confirmed by direct 
optical and infra red HST images (jCourbin et all l2002t 
iWinn et al.l[200l . 

A key species in the physics and chemistry of the 
ISM in galaxies is atomic carbon. For instance, [C I] 
492 GHz is one of the major far-infrared (FIR) fine struc- 
ture atomic lines via which cooling of the gas occurs 
(jHoUenbach fc Tielend [19971) . Moreover, both C and 
C"*", which directly follow from C, are essential in radi- 
cal rea ction networks, especially those leading to c arbon 
chains (jSuzuki et al.lll992t [Stonberg fc Dalgarnairi995) . 
Atomic carbon is also crucial for better understanding 
the structure of photo-dissociation re gions (PDRs). In- 
deed, early observations of [C I] fe.g.. lKeene et allll985f) 
revealed more intense emission than originally predicted 
by the models. Observations and theory now seem to 
have been reconciled by co nsidering the effects of the 
clumpy structure of clouds (jHoUenbach fc Tielend 1 1997I . 
and references therein). 

Several studies have investigated the use of the neutral 
atomic carbon (C) to CO abundance ratio {Nq/Nco) to 
investigate the physical conditions of the gas. A num- 
ber of such analyses make use of UV absorption data, 
and hence are limited to diffuse clouds in our Galaxy. 
In order to measure Nq and A^co in a dense environ- 
ment, many studies have targeted SubmiUimeter emis- 
sion of these species. However, the lowest C transition 
lies at a rather high frequency {^Pi Pq at ~492 GHz), 
and so its detection requires the use of ground-based fa- 
cilities with sensitive receivers and/or a large collecting 
area, or the observation of a redshifted source. More- 
over, in these SubmiUimeter studies, the emission line 
data are averaged over large areas and suffer from radia- 
tive transfer and excitation effects (with a bias toward 
the densest regions). The eSMA ("extended" Submil- 
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limeter Array, iBottinelli et al.ir2008f) is particularly well 
suited for this kind of observation as it provides the high 
sensitivity and high spatial resolution required to observe 
the weak [C I] transition in PKS 1830. In this work, we 
present the first data obtained with the eSMA, which 
probe Nc/Nco directly for the first time in absorption 
in the dense molecular regions of an external galaxy. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. eSMA and SMA data 

The eSMA consists of the SMA^ array of eight 6-m 
antennas augmented by the nearby single dishes of the 
JCMTi" (15 m) and the CSO (10.4m). The collecting 
area of the eSMA is twice that of the SMA alone, 
providing increased sensitivity, in particular for the 
longest baselines, which include the JCMT and CSO. 
Moreover, the longest baseline of the eSMA is increased 
by 50% compared to the SMA alone, which can bring 
the angular resolution down to below O'.'S at 230 GHz, 
as reported here, and to below 0'.'2 at 345 GHz, which 
will be the preferred frequency for eSMA observations. 

Observations of PKS 1830 (aj2ooo = 18'^33'"39!889, 
^J20oo = -21°03'39'.'77) were carried out with the eSMA 
for seven hours on 2008 April 14 (with eight SMA anten- 
nas in the very-extended configuration, or "vex" ) , and for 
eleven hours on 2008 August 5 with the SMA alone, also 
in vex configuration. The eSMA observations targeted 
the C I (^Pi-^Po) transition while the SMA provided 
CO, i^CO and C^^O (all J = 4-3) data. Table[T]summa- 
rizes the main observational parameters. The 45 eSMA 
baselines and the 28 SMA baselines had lengths rang- 
ing from 25-782 m (22-678 kA) and 68-509 m (56-419 kA) 
respectively, resulting in a best resolution at 1.1 mm of 
0'.'55x0'.'22. The tracks interleaved 15 to 20- minute ob- 
servations of PKS 1830 with 3 to 5-minute observations 
of the calibrators. Observations were carried out in sin- 
gle polarization mode with half-wave plates in the beam 
of the CSO and JCMT that automatically rotate to a 
common polarization as a function of elevation. 

In this interim state of eSMA commissioning, the 
JCMT used receiver A3 which, due to a different in- 
termediate frequency, provided a usable bandwidth of 
about 1.5 GHz, slightly less than the 2 GHz covered by 
the SMA and CSO. For both eSMA and SMA tracks, the 
correlator was setup to provide a uniform spectral resolu- 
tion of 0.8125 MHz per channel across the full 1.5-2 GHz 
bandwith, resulting in 2304 to 3072 channels per base- 
line and velocity coverages of ~1725 and ^2450 km s~^, 
for the eSMA and SMA respectively. Weather conditions 
were poor on the two nights with atmospheric opacity at 
225 GHz ranging from 0.15 to 0.25. 

2.2. Data reduction 
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tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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Astronomy Centre on behalf of the Science and Technology Facil- 
ities Council of the United Kingdom, the Netherlands Organisa- 
tion for Scientific Research, and the National Research Council of 
Canada. 



The data were reduced and calibrated using the 
MIR/IDL software package^^, and imaging was carried 
out with the MIRIAD software package. The SMA re- 
ceiver temperatures ranged from ~45 to 70 K on both 
nights, with the CSO comparable to the high end of 
this range, and the JCMT a factor of 2.5 times higher. 
The IF passband phase and amplitude were calibrated 
with observations of the strong sources 3C273, 3C454.3, 
J1924-292, and NRAO530 for the eSMA observations, 
and 3C273, 3C279 and 3C454.3 for the SMA-only data. 
Self-cahbration of the eSMA phase of PKS 1830 was 
done on the shortest possible timescales using an ini- 
tial model consisting of a pair of equal point sources 
separated by 1'.'05 at a position angle of 42°. The am- 
plitude gains for PKS 1830 were determined from the 
self-calibrated amplitudes of J1924-292. For the SMA- 
only data, phase calibration was performed in a standard 
way using J1733— 130 to derive time-dependent ampli- 
tude gains throughout the night. The amplitude gains 
imply average efficiencies of about 0.75 for the SMA an- 
tennas, 0.70 for the JCMT, and 0.30 for the CSO; the low 
efficiency of the CSO is likely due to a known malfunction 
of its rotating polarizer. To set the flux scales, standard 
SMA monitoring observations of J1924-292 (eSMA) and 
3C279 (SMA-only) at 230 GHz and 345 GHz were used 
(uncertainty < 20%). 

3. RESULTS 

The continuum map, obtained with the eSMA and de- 
rived from line-free channels, is displayed in Fig.[TJ This 
map shows the two lensed images of the distant quasar, 
separated by 1", consistent with previously reported val- 
ues. Figure [T] also shows the [C I] spectra at the peak of 
the continuum emission of the SW and NE components. 

Figure [5] shows the absorption spectra of [C I] and 
CO (4 — 3) as a function of heliocentric velocity, Vhci, 
observed in the SW source by the eSMA and SMA re- 
spectively. Two- and three-component gaussian fits were 
apphed to the [C I] and CO spectra respectively and are 
overplotted on Fig. [2l Our [C I] velocities and widths 
agree within the er rors with tho s e rep o rted for other 
mole cular species (jMuller et al.l 120061 : iMenten et al.l 
l2008f l and these values were used to guide the fits 
to the blended CO components. The presence of a 
fourth velocity component at ^^+20 km s^^ in the CO 
spectrum is possible, but was not fitted due to its low 
signal-to-noise ratio and is therefore not investigated 
further in this work. The derived parameters are given 
in columns 2-5 of Table [H The optical depths at 
velocities T4ici,o were derived using tq = — ln(Fv4^j „/ F^), 
where -FVhd o -^c are the line and continuum fluxes 
respectively. This equation implicitly assumes that the 
covering factor /c of the absorbing molecular gas across 
the finite extent of the background continuum source is 
1, which will be discussed in Section [H 

Using line data from LAMDA, the Leiden Atomic 
and Molecular Databa se^^, we ran the radiative trans- 
fer code RADEX (van der Tak et an [2007) for a range of 
H2 number densities n and kinetic temperatures Tkin, 
and for Tcmb=5.16 K; we obtain the excitation tem- 
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perature for each (n, Tkin), and hence determine (n, 
Tki„)-(0.2 - 5 X lO^cm-s, 10 - 30 K) corresponding 
to the excitation temperatur e Tex=rrot'^5.2-8 K d erived 
by iMuller et all ^06) and iWiklind fc Combes! (|T996h 
for molecular diagnostics of the physical conditions of 
dense cores such as HCO+. Using these same (n, Tkin), 
we in turn obtain Tex for CO and C of ^^6-20 and 8- 
25 K. These slightly higher Tox values are consistent 
with the lower Einstein A-coefficients of CO and C com- 
pared to HCO+. A ^co and Ng ar e then computed from, 
e.g., Eq. 2 of Mu Uer et all (|2006D . and are given in Ta- 
ble [2] The largest Nc is co nsistent with the up per limit 
of ^ 10^^ cm~2 reported bv lGerin et al.l (jlQQTn . and the 
CO column density we obtain for the ~7 km com- 
ponent, 0.9^0 5 X 10^^ cm~^, is also in agreement with 
the estimate given by these authors. Our results yield 
A^c/A^co'^2.5 for the two components at T4ici,o ^ —7 
and ^6 km s^^, and an upper limit of ^0.5 for the third 
component at Vtici,o — 19 km s^^. 

Finally, we derived an upper limit on A^co in the NE 
component of 5, 3.5 x 10^*^ cm~ ^ , assuming a line width of 
15 km s~i ()Muller et al.ll2006f) . This is consistent with 
the smallest A^co of 3 x 10^^ cm~^ expected from the 
range of HCO"*" to CO colur nn density rat i os det ermined 
from our data and that of iMuller etall (|2006f) . Non- 
detection of the targeted molecules in the NE component 
could be due to the ch ange in absorption reported by 
IMuller fc GuelinI (|2008f) in this source. 

4. DISCUSSION AND CONCLUSION 

Observations of C, C+ and CO can be used to probe 
the physical conditions of the gas. However, studies of 
the (clumpy) interstellar medium in external galaxies 
are complicated by several issues: (i) the large observ- 
ing beams which encompass several regions of differing 
conditions; (ii) the difhculty of performing [C I] obser- 
vations due to the weakness of the lines and the poorer 
atmospheric transparency at the high frequencies of the 
[C I] transitions (492 and 809 GHz); and (iii) the impos- 
sibility to access the [C II] 158-/im transition from the 
ground. The observations presented here make it possi- 
ble to tackle the first two points. Indeed, since the lines 
are detected in absorption and the galaxy responsible for 
the latter is almost face-on, the line profiles will not be 
affected much by rotation broadening. Regarding point 
(ii), it is worthy to note that in fact, no [C I] absorp- 
tion has yet been detected in a local galaxy at z ^ 0. 
For PKS 1830, the absorbing material is at a redshift 
of -0.89, shifting the 492-GHz [C I] hne to ~261 GHz 
where the atmosphere is more transparent and stable. 

While absorption data circumvent some of the prob- 
lems inherent to emission spectra, some limitations un- 
fortunately remain, in particular, our (lack of) knowledge 
of the structure of the absorbing gas. This translates 
partly into the estimation of the covering factor fc, which 
is one of the major sources of uncertainty (along with 
Tex) in the determination of the column densities. For 
C0(4— 3), the largest optical depth is ~2, meaning that 
the absorption almost reaches the zero intensity, but it is 
not heavily saturated as indicate d by the non-detectio n of 
the less abundant ^^CO. Indeed. [Mrdler et all (l2006l ) de- 
rive ^^C/^^C ~27, so that if CO were heavily saturated, 
the expected optical depth of ^'^CO would be >0.1, i.e. 
at the 2-cr detection limit. Hence, the covering factor of 



the absorbing material must be close to unity, as noted 
by previous studies of this source. This is consistent with 
the picture of indi vidual clouds with sizes of order ~1 pc 
(cf. (Menten et a"Llll999l : IMuller fc Gu6 1in 2008, and sizes 
of typical Galactic molecular clouds of up to a few pc) 
in front of a background source of projected size <1 pc 
at the distance of the galaxy, corresponding to the de- 
convolved core size of ~ 0.228 X 0.148 mas measured by 
iJin et all (|2003l ) for the SW image of PKS 1830 at 7 mm 
with the VLBA. 

Comparing the total H column density derived with 
different met hods can indicate whether the "cloud" is 
uniform (e.g. IWiklind fc Combes! [l99l . For the SW 
source, there is no strong evidence for dumpiness, since 
in this case the total H column density derived from 
X-ray data (- 3.5 x 10^^ cm^^ Mathur fc Nair 1997) is 
comparable to that derived from millimeter CO absorp- 
tion ('^ 5 X 10^^ cm~^, assuming C0/H2~ 5 x 10~^), 
or from C I absorption ( ^ 6 x 10^^ cm~^, as suming a 
typical C/H2~ 8 X 10-^ iFrerking et 31111989( 1. Hence, 
the assumption that the individual velocity components 
consist of material with homogeneous physical and 
chemical conditions seems sufficient, which allows us to 
compare the column density ratios. 

Our data indicate two types of Nc/Nqo (see Table 3 
for comparisons): a low value, ^0.5, representative of 
dense cores or PDRs (< 0.5 — 1), and a somewhat high 
value, ~2.5, close to those derived for translucent clouds 
(~ 3 — 6). Qualitatively, low Nq/Nqo are found in en- 
vironments with low neutral atomic carbon abundances, 
such as high-UV environments (where C becomes ion- 
ized), or high-density environments (where CO forma- 
tion is efhcient). The presence of high density tracers 
such as HCO+, HCN and CS supports the latter case 
of a dense core for the Vhoi = —19 km s~^ component. 
On the other hand, high Nq/Nco are representative of 
low-density, low-column density, mild-UV environments 
where CO is photodissociated but atomic carbon is not 
ionized. For the Vhd ~ —7 km compone nt, the de- 
tection of hot NHg (VisR km s"\, Hcnkcl eFaL|200l) 
and l arge H I optical depth (e.g. iKoopmans fc de Bruvnl 
l2005l ) point toward a low density PDR. Since hot NH3 
was not detected at 14ei ~ 5 km s~^, the high Nc/Nco 
here more likely indicates diffuse or translucent condi- 
tions. An absorbing source at z ~ 2.4 with similar diffuse 
m aterial, albeit a higher Nq/Nqo of ~10, was reported 
bv lSrianand et al.l (2008. ). 

Atomic C variati ons could also be e xplained by 
chemical evolution (jMaezawa "eTall [19991) : since the 
timescale for conversion of C to CO is comparable to 
the dynamical timescale of a dense core 10^ yr), 
Nq/Nco is expected to be high in the early stage of 
dense core formation and to decrease with time. Hence, 
our 7Vc/A"co~2.5 could be indicative of a cloud in an 
intermediate stage of evolution from diffuse to dense 
gas, where C has not completely converted into CO. 

In conclusion, the data presented here demonstrate the 
ability of the eSMA to provide observers with high an- 
gular resolution, high sensitivity observations. Further- 
more, they allowed us to investigate the physical and 
chemical conditions of the material obscuring the SW 
component of PKS 1830, and determine Nc/Nqo di- 
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rectly in a dense molecular cloud located at ^=0.886. 
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TABLE 1 

Main observational parameters for PKS 1830 



Parameter 


DSB continuum 


C ePi-^Po) 


i2CO(4-3) 


I/rost (GHz)'' 




492.16065 


461.04077 


!^obB (GHz)=' 


1.1 mm 


260.97965 


244.47761 


Channel width 


2xr-1.5 GHz'^ 


0.93 km s-i<: 


0.99 km s-i<: 


Beam FWHM 


0'.'55 X 0('22 


0'.'62 X 0'.'30 


0'.'57 X 0('40 


Beam P.A. 


32° 


38° 


15° 


Array 


eSMA 


eSMA 


SMA 



TABLE 2 

Observed and derived line parameters for the SW component 



Transition 


Vhol.O 

(km s~^) 


To 


FWHM 
(km s"^) 


frdV 
(km s"^) 


AT" 

(lO" cm^^) 


JVc 
A^co 


[C I] 


-19.0 


<0.15= 


10. o" 


<1.6 


<1.6 


<0.6 




-7.8 


0.65 


15.1 


10.3 


10.21:^1 


^- ' -2.0 




5.0 


0.99 


12.9 


13.6 




n O + 0.4 
^■^-1.6 


CO 


-19.0*' 


0.86 


13.7 


12.5 


1.3 






-6.5'' 


1.59 
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18.6 


1.9 
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2.21 


13.0'' 


30.6 
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"CO, c^**o 




<0.1'' 


10. 0"* 


<1.1 
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TABLE 3 

Comparison of Nq/Nco measurements. 



Source 


Nc/Nco 


References 


Translucent clouds 


~3-6 


Stark &; van Disliocck 1994 


Diffuse clouds 


~l-80 


Fcderman et al. 1980 


Diffuse material at 2 — 2.418 


~10 


Srianand et al. 2008 




Molecular 


clouds 


Ophiucus 


~0.05-l 


Frcrkine et al. 1989 


G34.3-I-0.2 


~0.1-1 


Little et al. 1994 


W3, NCG2024, S140, Cep A 


--0. 1-0.5 


Plume et al. 1999 


TMCl, L1527 


~0. 06-0. 25 


Maezawa et al. 1999 


Orion A/B (multi position) 


0.05-0.21 


Ikeda et al. 2002 


PDRs (multi positions) 


M17, S140; IC 348 


0.02-0.25 


Keene et al. 1985: Sun et al. 2008 


NGC1977 


0.03-0.2 


Minchin & White 1995 


Orion Bar 


0.04-0.3 


Jansen et al. 1995 


S106 


0.02-0.09 


Schneider et al. 2003 


IC 63 (high/low density) 


0.1-0.7 / 2-3 


Jansen et al. 1996 


Galaxies 


M33 (4 positions) 


0.03-0.16 


Wilson 1997 


IC 342 (center) 


0.5-0.6 


Israel & Baas 2003 


NGC 253; Maffei 2 (center) 


0.2-0.3 


Israel et al. 1995; Israel & Baas 2003 


M82 


~0.5 


White et al. 1994 


M83 (5", 5") 


~0.3 


Petitpas & Wilson 1998 


M83 (center), NGC 6946 


~0.9 


Israel & Baas 2001 


M51 (center) 


0.4-1.0 


Israel et al. 2006 


Spiral at 2=0.886 


<0.5, ~2.5 


This work 
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Bottinelli et al. 
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Fig. 1. — eSMA 1.1 mm continuum map and [C I] spectra of PKS 1830. Absolute astrometry information is lost during the self-calibration 
process, so that the offsets here are relative to the center of the double point source. The contour levels are multiples of 150 mjy. The 
0'.'55x0'.'22 beam is shown in the lower left corner. 
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Fig. 2. — eSMA spectrum of C I (^Pi — ^ Po) and SMA spectrum of CO ( J = 4 — 3) observed towards the SW component. Spectra were 
smoothed to a velocity resolution of 2.8 and 5.9 km s~^, yielding rms before normalization of 0.07 and 0.04 Jy beam~^, corresponding to 
0.09 and 0.05 on optical depth scale, for [C I] and CO respectively. Overlaid in dashed lines are two- and three-component gaussian fits to 
the [C I] and CO profiles respectively, while the solid lines are the sum of these individual components. 



